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Abstract 

We present a detailed study of the structural, electronic, magnetic and ferroelectric properties 
of prototypical proper and improper multiferroic (MF) systems such as BiFe03 and orthorhom- 
bic HoMnC>3, respectively, within density functional theory (DFT) and using the Heyd-Scuseria- 
Ernzerhof hybrid functional (HSE). By comparing our results with available experimental data as 
well as with state-of-the-art GW calculations, we show that the HSE formalism is able to account 



c/5 



O . 

well for the relevant properties of these compounds and it emerges as an accurate tool for pre 



' dictive first-principles investigations on multiferroic systems. We show that effects beyond local 
■ and semilocal DFT approaches (as provided by HSE) are necessary for a realistic description of 

OS 

MFs. For the electric polarization, a decrease is found for MFs with magnetically-induced ferro- 

m ; 

electricity, such as HoMnOs, where the calculated polarization changes from ~ 6 fiC/cm using 
Perdew-Burke-Ernzerhof (PBE) to ~ 2 fiC/cm 2 using HSE. However, for proper MFs, such as 

o ■ 

BiFeC>3, the polarization slightly increases upon introduction of exact exchange. Our findings 
therefore suggest that a general trend for the HSE correction to bare density functional cannot be 
extracted; rather, a specific investigation has to be carried out on each compound. 



PACS numbers: 75.80.+q; 75.47.Lx; 75.50.Ee; 77.80.-e; 71.15.Ap; 



I. INTRODUCTION 



Multiferroics (MFs) are materials in which different ferroic orders such as ferromagnetism, 

n 

ferroelectricity and/or ferroelasticity may coexist in a single compound. [1] They have at- 
tracted much attention for their potential applications in memory devices and other elec- 
tronic components, due to the intriguing possibility of controlling magnetism by an applied 
electric field, and viceversa (magnetoelectric effect). {2 5 1 

Multiferroics are compounds where electron correlations are rather important, and where 
the electron charge shows atomic-like features, such as strong space localization, poorly 
dispersed band energies, and large on-site Coulomb repulsion. 

aaaa 

For these sys- 
tems, there are well-known deficiencies of local-spin-density approximation (LSDA) or spin- 
polarized generalized-gradient-approximation (SGGA) to density-functional-theory (DFT). 
Among them, we recall the underestimation of the band-gap magnitude for most insulating 
materials, jg- 10 1 Part of these failures can be traced back to the self- interaction error in ap- 



proximate density functionals: the electron charge experiences a spurious interaction with 
the Coulomb and exchange-correlation potential generated by itself. 
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The LSDA+C/ 
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13Ml5| and the self-interaction correction (SIC) schemes 



16 



- |l8| can over- 



come some of the deficiencies of LSDA. However, LSDA+[7 suffers ambiguities in the choice 
of the U parameter and needs a choice regarding which orbitals to treat within a Hubbard- 
like approach. For simple materials, a self-consistent evaluation of the U parameter can be 
obtained, although this method is not widely used.[19( For BiFe0 3 , the value U e //=3.8 eV 



has been recently calculated. 



20| 



SIC-schemes are not commonly available in electronic structure codes for extended 



solid state systems. The implementation of a fully sel 



-consistent SIC-LSDA approach 



cers. 2lj Since then, other approaches 



22]). SIC-schemes suffer from the 



for extended systems was done by Svane and co-worl 
have been implemented (for a review, see Refs. 
"nonvariationality-problem" of the energy functional which makes forces and stress cal- 
culation not commonly available. 



23] 



In the last few years, h ybr id Hartree-Fock density functionals 24j-l27f have 



3een widely 



used in solid state physics. [27H30I] ranging from simple semiconductor 



sition metals, lanthanides, actinides. 
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341 molecules at surfaces. 
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systems, 



3CH ] to tran- 



semiconductors,|37] carbon nanostructures. 



38 



36| diluted magnetic 



For a recent review see Ref. 



40|. Hybrid 



functionals mix the exact nonlocal exchange of Hartree-Fock theory 24j-l27j with the density 
functional exchange. The Heyd-Scuseria-Ernzerhof screened hybrid functional (HSE) 4ll. |42| 
is well suited for extended solid state systems. 

There are very few studies dealing with ferroelectric oxides and even less with multi- 
ferroics. Wahl et al. re-investigated the well-known SrTiC>3 and BaTiC>3 [43y using HSE 
and semilocal functionals (LDA,PBE,|44] PBEsol|96]). Bile et al. studied in great details 
BaTiC>3 as well as PbTi03 using the Bl-WC hybrid-functional and concluded that the latter 



gives an accurate description of both the structural and electronic properties. |45J Goffinet 
et al. extended the analysis to the prototypical multiferroic bismuth ferrite showing that 
hybrid-functionals, specifically the Bl-WC functional, open new perspectives for a better 
first-principles description of multiferroics. 46( In passing, we recall that the WC gradient 
corrected functional [47] is very similar to PBEsol and the hybrid Bl-WC 48| functional mixes 
the WC functional with 16% nonlocal exchange. The hybrid HSE functional mixes 25% non- 
local exchange with the PBE functional and the mix is performed only on the short range 
component of the Coulomb interaction (for further details, we refer to Ref. |43T ). However, 
which functional to prefer for simple ferroelectric compounds is still an open issue. [43] 

So far, a good performance of HSE or Bl-WC functionals has been recognized for proper 
multiferroics where the ferroelectric polarization is of displacive type. On the other hand, 
magnetically driven multiferroics, also known as improper multiferroics, are largely un- 
explored using hybrid functionals. The purpose of this work is to extend the previous 
hybrid density functional studies from prototypical ferroelectric oxides (SrTi0 3 , BaTi0 3 , 
PbTiOs) 43|, |45] or simple multiferroic system (BiFeOs) to more complicated and exotic 
multiferroic compounds, such as HoMn0 3 . 

First of all, we focus on BiFe0 3 , already investigated using the Bl-WC functional, but 
not yet using HSE. In this way, we are able to compare two different, although similar, 
approaches for BiFe03. Most importantly, we consider another prototypical case of improper 



multiferroic, namely HoMnOs, which has recently attracted much attention. |49l. |50| We will 
show that important differences compared to standard DFT approaches arise when a proper 
description of correlated electrons, such as that given by HSE, is taken into account. 

Our study suggests that HSE functional improves the description compared to standard 
DFT approaches for multiferroic systems. 

The paper is organized as follows. Details of the computational setups are given in 



Sect. [TTl An extended discussion of the structural, electronic, magnetic and ferroelectric 
properties of BiFeC>3 is reported in Sect. IIHI Sect. [IV] is devoted to HoMnC>3 focussing on 
the paraelectric AFM-A (Sect. HVAj) and ferroelectric AFM-E (Sect. HVBj) ph ases. Finally, 



in Sect. |Vl we draw our conclusions. 



II. COMPUTATIONAL DETAILS 



All the calculations presented in this study are performed by using the latest version of the 
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For BiFe0 3 , all the results are obtained 
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53| by explicitly treating 15 valence 



Vienna ah initio simulation package (VASP 5.2). 
using the projector-augmented plane-wave method 
electrons for Bi (5<i 10 6s 2 6p 3 ), 14 for Fe (3p 6 3<i 6 4s 2 ), and 6 for oxygen (2s 2 2p 4 ). We used a 
6x6x6 Monkhorst-Pack k-mesh for the Brillouin-zone integration and 400-eV energy cutoff. 
Tests using a 8x8x8 mesh as well as 600 eV cutoff did not give significant differences in the 
calculated properties. Brillouin zone integrations are performed with a Gaussian broadening 
of 0.1 eV during all relaxations. The experimental unit cell for the R3c (ferroelectric phase) 
was used as an input in the full-optimization procedure. For this phase as well as for 
the paraelectric one (see below), we used the rhombohedral setting. The geometries were 
relaxed until all force components were less than 0.01 eV/A and the stress tensor components 
less than 50 meV/cell. The spin configuration was fixed in order to reproduce the G- 
type antiferromagnetic state of BiFeOa and the spin-orbit coupling was neglected. For the 
paraelectric phase, we used the non-polar R3c LiNbOs phase. [54] We compute the difference 
of electric polarization, i.e. AP = P FE - P PE = (PfJ + P™) - (PfJ + P™) = AP ton + 
AP e i e , where FE, PE, ion and ele denote ferroelectric, paraelectric, ionic and electronic 
contribution, respectively. For the paraelectric phase, we used the same lattice constant and 



PE 



rhombohedral angle of the ferroelectric one. Note that, although counterintuitively, P e ^ e 



may be different from zero, as explained in Ref. [55|. Pf^ PE is calculated by summing the 



position of each ion in the unit cell times the number of its valence electrons. The electronic 



contribution is ob 
of polarization. 56 



gained by using the Berry phase formalism, within the "modern" theory 



Concerning the HSE calculations, due to the high computational load,[97| we always used 
the 400 eV and 6x6x6 k-point mesh. The Fock exchange was sampled using the twofold 
reduced k-point grid (using the full grid, gives however negligible changes in the computed 



properties). Finally, we performed GqWq calculations 59H64j on top of the HSE electronic 



and ionic structure, which usually represent a good starting point for a perturbative quasi- 



particle excitation energies. 



60| We also included vertex correction in W via an effective 



nonlocal exchange correlation kernel. 63| 



For orthorhombic HoMnC>3, the Pnma symmetry is chosen with the b basis vector as the 
longest one. The paraelectric phase was simulated using 20-atoms cell in the AFM-A spin 
configuration showing ferromagnetic (FM ) (AFM) intraplanar (interplanar) coupling; for the 
ferroelectric one we used a 40-atoms cell (doubling the previous cell along the a axis) in the 
AFM-E spin configuration (i.e. in-plane FM zigzag chains anti-ferromagnetically coupled to 
the neighboring chains with the interplanar coupling also AFM). The energy cutoff was set 
to 300 eV and the Brillouin zone mesh was fixed to 4x2x4 and 2x2x4 grid for the AFM-A 
and AFM-E phase respectively. Ho 4/ electrons were assumed as frozen in the core. The 
experimental lattice constants were used for all the calculations but the internal positions 
were relaxed. For the HSE calculations, the Fock operator was evaluated on the down-folded 
k-point mesh. In order to assess the relative stability of the two magnetic phases, we used 
the same simulation cell containing 40 atoms for both phases, increasing the cutoff to 400 
eV and using a 4x2x4 k-point grid. 

III. BIFEO3: RESULTS AND DISCUSSIONS 
A. Structural properties 

The ferroelectric structure is represented by a distorted double perovskite structure with 



R3c symmetry (N. 161, point group C3„) as reported by Kubel and Schmid. 65[ The para- 
electric phase has R3c symmetry (N. 167, point group D 3 d). Both phases are shown in FigJH 
In Table [I] we report relevant properties such as the structural parameters, the Fe magnetic 
moment and the energy gap calculated using the PBE and HSE functionals. We also report 



the values using the Bl-WC functional taken from Ref. 46] 

First of all, HSE reduces the lattice parameter a r h with respect to PBE, giving a much 
better agreement with the experimental value: the error decreases from ~1 %(PBE) to 
~0.3 %(HSE). As a consequence, the unit cell volume V also shrinks, getting closer to 
the experimental value. The rhombohedral angle, a r h, is almost insensitive to the applied 



functional. Thus, the inclusion of Fock exchange makes the structure more compact, i.e. 
the lattice constant decreases. Note that the Bl-WC functional gives too small lattice 
constant and too small equilibrium volume as compared to HSE, worsening the comparison 
with the experiments. There is a very good agreement between the relaxed coordinates of 
the Wyckoff positions and the experimental ones using HSE, while the PBE as well as the 
Bl-WC functional give slightly worse results ( the only exception being the x component 
of the oxygen atoms in the 6b site symmetry). In the experimental structure, the BiOe 
cage is strongly distorted with three coplanar nearest neighbors (NNs) lying above Bi along 
[111] at 2.270 A (d' m _ , s refers to short) and three NNs sitting below at 2.509 A (d l Bi _ , I 
refers to long). From Table HI we see that the theoretical NNs distances compare well with 
experiments, with errors from ~ 1 % to ~4 % (PBE), from ~ 3 % to ~ 5 % (Bl-WC) and 
~ 2 % (HSE). The O-Fe-0 bond angle would be 180° in the ideal cubic perovskite. In this 
system, it buckles to an experimental value of 165.03°. The HSE value (164.56°) is close to 
PBE, and in both cases, they are slightly underestimated with respect to experiment. The 
Bl-WC angle, on the other hand, is clearly underestimated. Overall, the predicted HSE 
values clearly are in much better agreement with experiments than those calculated using 
the PBE or Bl-WC functional. 

B. Electronic and magnetic properties 

Let us consider now the magnetic and electronic properties. As shown in Table HJ the 
calculated local moments are generally very similar for all the functionals, and close to 
the experimental value. In particular, the HSE local moment is slightly larger than PBE 
, suggesting a more localized picture of the spin-polarized electrons. The calculated PBE 
(HSE) electronic energy gap is 1.0 (3.4) eV. The expected band-gap opening using hybrid 
functionals can be understood as follows: the exact exchange acts on occupied states only, 
correcting them for the self-interaction, thus shifting downwards the occupied valence bands. 
In turn, this has a clear interpretation: within the Hartree-Fock approximation for the 
ground state of an N electron system, the potential felt by each of the N electrons in the 
ground state is that due to N-l other electrons, i. e. they feel a more attractive ionic potential. 
On the other hand, for unoccupied states, the potential is that due to the N occupied orbitals, 
so these orbitals effectively experience a potential from one more electron, the latter screens 



the ionic potential which in turn becomes less attractive. Therefore, the unoccupied states 
are shifted upwards, opening the gap. 

As for the experimental energy gap for BiFeOs, the situation is not clear. There have 
been several measurements of the band gap using UV-visible absorption spectroscopy and 
ellipsometry on polycrystalline BFO films, epitaxial BFO films grown by pulsed-laser depo- 
sition, nanowires, nanotubes, and bulk single crystals. Reported band-gap values vary from 



2.5 to 2.8 eV. 
Kanai et al 
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7l| An estimate gives ~ 2.5 eV from the optical absorption spectra by 
70j and Gao et al. [ri]. From the theoretical side, there is a spread of values: a 



small gap of 0.30-0.77 eV using LSDA,[55j or from 0.3 to 1.9 eV using "LDA+U" , depending 



on the value of U 



55|; 0.8-1.0 eV usingPBE (WC) GGA functional; |45| 3.0-3.6 eV using Bl- 



WC and B3LYP hybrid functionals. 45| Thus, a parameter-free theoretical reference value 
is clearly needed. The most accurate (but expensive) method is the GW approximation. js9 ] 
Here, we provide for the first time, the value of the BFO energy gap based on the GW 
method. 

First of all, at the PBE level, we estimate an energy gap of 1.0 eV. When introducing the 
exact exchange (HSE), the gap opens up to E g =3.4 eV. Upon inclusion of many-body effects 
(GoWo), it opens even more (E g =3.8 eV). Finally, when including vertex corrections, we find 
that the gap reduces to 3.3 eV. Remarkably, vertex corrections almost confirm the HSE band 
gap. This is perfectly in line with recent works 40j where it is argued that HSE band gaps 
represent a very accurate estimate due to partial inclusion of the derivative discontinuity 
of the exchange-correlation functional. [8] Clearly, effects beyond bare DFT are important 
in this compound. Our results show that hybrid-functional calculations give already a very 
good estimate at a lower computational cost compared to GW. In this respect, we mention 



72| 



the very recent experimental study based on resonant soft x-ray emission spectroscopy 
where the band gap corresponding to the energy separation between the top of the O 2p 
valence band and the bottom of the Fe d conduction band is 1.3 eV. The discrepancy between 
theory and experiment may be due to the presence of defects in the experimental sample as 
well as to the resolution involved in photoemission spectra. We hope to stimulate further 
experimental work to test our first-principles prediction of the energy gap of this important 
multiferroic material. 

In Fig. [2] we show the Density of States (DOS) for the optimized atomic structure. Let's 
focus on the PBE DOS. The lowest states at ~ — 10 eV are Bi s hybridized with O p states 



(blue curve). Above — 6 eV there are hybridized O p and Fe d states. The Fe d states extend 
in the conduction band as well; the Bi p states can be found above 4 eV. As for the HSE 
DOS, we see that the conduction bands are shifted upwards, opening the valence- conduct ion 
gap. There is a change in the spectral distribution above — 8 eV: a valley appears around 
— 6 eV and the lower shoulder of the peak increases its spectral weight. It is easy to trace 
back the above changes to modifications of the majority Fe d states, as shown in the panel 
beneath: while in PBE the band states in the vicinity of the top of the valence band have 
predominantly d character, in HSE the spectral weight of the Fe d states is concentrated far 
away from the top of the valence band. This can be interpreted as a change from a more 
itinerant picture to a more localized description of the Fe d states going from PBE to HSE. 
In Fig. |2l we also include the spectral distribution of the Fe d states derived from a recent 
experimental work 72j (see dotted lines): the position of the main HSE Fe d peak almost 
perfectly matches the experimental PDOS, although the bandwidth of the calculated DOS 
is different because of energy resolution, etc. Indeed, if we include the Fe d DOS calculated 
using G0W0+ vertex corrections, the agreement between the theoretical and experimental 
peak position becomes excellent. Note that the HSE and Go Wo + Vertex peak position are 
very close to each other, confirming the accurate HSE description of the BiFeOs electronic 
structure. In passing we note that while DFT+C/ gives a better estimate of the Fe d peak 
position, [55] the energy gap is still underestimated with respect our GW calculation. 

We previously mentioned that HSE may change the ionic/covalent character in this com- 
pound. To support this, in Fig. [3] we show the difference between the PBE and HSE charge 
density, Ap = p PBE — p H calculated at fixed geometry. In a purely ionic description, all 
the valence electrons would be located on the oxygens, acting as "electron sinks", and the 
cations would donate their nominal valence charge. The more the electrons populate the 
anions, or conversely, the more the electrons depopulate the cations, the more the picture 
shows an ionic character. Fig. [3] confirms the trend discussed before: upon adding a fraction 
of exact exchange to the PBE functional, the electronic charge at the cations decreases, i.e. 
Ap is positive (grey areas in Fig. [3]), while the electronic charge at the anions increases i.e. 
Ap is negative (yellow areas in Fig. EJ). Thus, the introduction of exact-exchange generates 
a flux of charge from the cations towards the anions, clearly shown in Fig. El increasing the 
ionicity of the compound. 

In order to discuss more quantitatively these effects, we perform a Bader analysis of 



the electronic charge. 73H75j The atom in molecules (AIM) theory is a well established 
analysis tool for studying the topology of the electron density and suitable for discussing 
the ionic/covalent character of a compound. The charge (Qb) enclosed within the Bader 
(Vb) volume is a good approximation to the total electronic charge of an atom. In Table 
Ull we report Q B and Vb calculated for Bi, Fe, and O at a fixed geometric structure, i.e. 
HSE geometry. This is needed in order to avoid different volumes for the normalization of 
the charge in the unit cell and for highlighting the electronic structure modifications due 
to the exact exchange. Furthermore, we consider only the valence charge for our analysis 
(although one should formally include also the core charge, we do not expect variations 
as far as the trends are concerned). Let us first consider the cations: the Bader charge 
and volume are larger in PBE than in HSE. For the anions, the opposite holds true. This 
is not unexpected and in agreement with intuition: upon introducing Fock exchange, the 
system evolves towards a more ionic picture, through a flux of charge from cations towards 
anions, which reduces (increases) the "size" of the cations (anions) when going from the 
PBE to the HSE solution. Finally, we note that a different degree of ionicity modifies the 
calculated equilibrium lattice parameter: in a partially covalent material, such as BiFeO^.f?^ 
the increased ionicity changes the different net charges generating a higher Madelung field, 
which is an important contribution to the bonding in the solid, and contracts the equilibrium 



structure. 
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C. Ferroelectric properties 

Let's finally focus on the electric polarization. In Table HVt we report the ionic and 
electronic contributions to the difference of ferroelectric polarization between the polar (R3c) 
and non-polar (R3c) both in PBE and HSE. In order to disentangle the purely electronic 
effects from the ionic ones upon introduction of Fock exchange, we report also the PBE(HSE) 
electronic contribution calculated at fixed HSE(PBE) geometry. 

As a general comment, we note that a large polarization of ~ 100 fiC/cm 2 along (111) 
for bulk BFO has been reported experimentally by new measurements on high-quality single 
crystals. |79| in good agreement with our calculated values. In what follows, we will mainly 
focus on the differences between PBE and HSE calculations. Note that the unit cell volume 
is different for PBE and HSE, as shown in Table [H 



First, we note that the polarization calculated according to the point charge model (P pcm ) 
is closer to ~P to t at the HSE than at the PBE level. This confirms that the HSE description of 
BFO points towards an ionic picture, i. e. by decreasing the covalency effects. The calculated 
total polarization P io4 is ~ 105 fiC/cm 2 using PBE and ~ 110 fiC/cm 2 using HSE, i.e. 
HSE predicts an increase of total electric polarization. The occurrence of ferroelectricity in 
BiFeOa is usually discussed in terms of "polarizable lone pair" carried by the Bismuth atom. 
This has a physical interpretation in terms of cross gap hybridization between occupied O 



2p states and unoccupied Bi 6p states. 



80h83|. Intuitively, the larger the energy gap, 



the lower the polarization should be. Accordingly, one might expect HSE to reduce the 
polarization compared to PBE because of the larger energy gap. We will show below that 
this is not in contraddiction with the results of Table [B In fact, let's consider P tot calculated 
at the same atomic structure (for example at the PBE relaxed structure of the paraelectric 
and ferroelectric phases) but using both PBE and HSE. We denote the former as Pf D f E , 
the latter as ¥^ e<kPBE \ Note that the ionic contributions is of course the same for both 
cases. From Tabled! we have P£f s =105.6 /iC/cm 2 and Pf a t E — 103.2 jjC/cm 2 , i.e. 
a decrease of total polarization is found in going from PBE to HSE for the same ionic 
structure. A similar behavior is found in opposite conditions: for the HSE ionic structure, 
P^f- E =110.3 /iC/cm 2 and P P ) f E ^ HSE ^ = 112.6 /jC/cm 2 . Thus, keeping the same volume and 
including Fock exchange, the polarization reduces as expected. On the other hand, when we 
evaluate the total polarization at the appropriate equilibrium and relaxed structures using 
PBE and HSE , the ionic contribution also varies and one loses a direct connection between 
the increase of the energy gap and the decrease of total polarization. In our case, the total 
polarization, when evaluated at the appropriate equilibrium volume for each functional, 
increases from PBE to HSE. This clearly points out a strong volume-dependence of the 
polarization, therefore calling for a correct estimate of the volume (as provided by HSE). 



IV. HOMN0 3 :RESULTS AND DISCUSSIONS 



A. Paraelectric AFM-A phase 

1. Structural properties 

An extended review of the main properties of orthorhombic HoMnC>3 within a standard 



PBE approach can be found in Ref. [84j where it is also shown that the inclusion of the U 
correction worsens the structural properties. Therefore, in this paper, we will focus on the 
comparison between the predictions of HSE with respect to PBE results. In Fig. H]we show 
the perspective view of HoMnC>3 and the paraelectric (AFM-A) and ferroelectric (AFM-E) 
spin configurations in the c — a plane. 

In Table II III we report the optimized structural parameters in the AFM-A magnetic 
configuration, calculated using PBE and HSE. The in-plane short (s) and long (/) Mn- 
oxygen bond-lengths get closer to experimental values using HSE; on the other hand, the 
out-of-plane length is slightly overestimated with respect to the PBE and the experimental 
value. In order to quantify structural distortions, the Jahn- Teller (JT) distortion vector 
Q — [Q\iQy\ — [Vl — s, y^|(2m — I — s)] is often introduced. From Table IHIl it is thus 
clear that the HSE functional improves the JT distortion upon the PBE description: the 
magnitude of Q is Q=0.55 A and 0.61 A using PBE and HSE, respectively, whereas the 
experimental value is 0.59 A. As far as the structural angles are concerned, we first notice 
that the GdFe03-like tilting (a) in the Mn-06 octahedron is slightly overestimated using the 
hybrid functional: the deviation from the experimental angle is 3.7 (PBE) and 5.0 % (HSE) 
with the Mn-O-Mn in-plane angles calculated using HSE slightly reduced with respect to 
PBE. We note, however, that the experimental uncertainty on the angles may be up to 



1°, 



85] due to synthesis problems of ortho-HoMnOs. 50|, |85) The octahedral tilting is 



related to the ionic size of the rare-earth ion: 
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87] the tilting increases when the radius 



of the rare-earth atom decreases (for example, from La to Lu in the manganites series) . [84 ] 
In this respect, the tendency towards a larger octahedral tilting, upon inclusion of exact- 
exchange-functional, goes hand in hand with the reduced ionic size of Ho ion when going 
from PBE to HSE. As in the previous case, we performed a Bader analysis of the valence 
charge distribution. Results are shown in Table HT1 as expected, the "size" of the Ho ion 
is reduced within HSE. According to our previous discussion, the ionic/covalent character 



of the charge density is modified by HSE in favor of a more ionic picture. This will have 
important consequences for the electronic polarization, as shown below. 



2. Electronic properties 

In Fig. [5] we show the band structure for the AFM-A phase as calculated using standard 
PBE (left panel) as well as HSE (right panel) along the main symmetry lines. The PBE 
band structure shows a small gap equal to ~ 0.2 eV. The bands below ~ —2 eV are mainly 
oxygen p states and those 2-3 eV below (above) the Fermi level are mainly spin-up (spin- 
down) Mn d states. There is also a considerable weight of the Mn d states in the oxygen 
bands near the top of the valence band. [3] 

The group of bands between — 1 and — 2 eV are mainly d xy ,d yz with some small weight 
of d z 2. The two bands just below the Fermi energy are mainly d^^-like with some d xz 
weight; at Y,S,Z they become degenerate. Higher in energy, between and 1 eV, there are 
two more bands showing a similar behavior, i.e. degenerate at Y,S,Z, almost degenerate 
along Z-R and with a similar overall band dispersion. Even higher in energy, the are the 
Mn minority states. From Fig. we can extract the JT splitting (Ajt of the "e 9 " states), 
the CF splitting (Aqf between "e g " and "^23"), and the exchange splitting (Aex between 
majority and minority spin states, evaluated at the S point, for simplicity). We thus have 
Aj T =1.05, A c ^=2.34, and A E x =l 2.75 eV. The comparison between the PBE and HSE band 
structure highlights some differences. First we note that the oxygen bands are slightly shifted 
to lower binding energy together with the Mn d xy ,d yz bands. On the other hand, the average 
position of occupied Mn d xz , d x 2_ y 2 bands remain almost unchanged, but the band-width 
increases. The latter effect is mainly shown by the lowest of the two bands, i.e. the d^-like 



band. As well known, 88J] Hartree-Fock hamiltonians naturally leads to larger band-width 
and down-shift of electronic states, even for simple homogeneous systems. Thus, the larger 
band-widths obtained by HSE can not be simply connected to a stronger hybridization, 
because it is an intrinsic feature of Fock exchange for every electronic state. Indeed, we will 
show below that the d — p hybridization is expected to decrease using HSE. 

Ajt is evaluated as 3.4 eV, larger than in the PBE case, suggesting a stronger local 
distortion related to the Jahn- Teller instability. The increase of the JT splitting is linked 
to the increase of the energy gap as well, which is now ~ 2.7 eV. The exchange splitting, 



A^x=4.6 eV, is also larger along with an increased local Mn moments with respect to PBE. 
Note, that the enhancement of Jahn- Teller distortion by HSE is not unexpected. In fact, 
it was previously suggested that the HSE functional is able to reveal the Jahn- Teller effect 
for Mn +4 through a symmetry broken solution giving rise to an orbitally ordered state and 



consequent Jahn- Teller distortion. 



37j | The larger A jj- is mainly driven by a purely electronic 



effect due to the inclusion of Fock-exchange. Infact, by calculating the PBE self-consistent 
charge density on top of the HSE ionic structure, Ajt becomes ~ 1.2 eV, i.e. nearly equal to 
the previous PBE case. Thus, we expect HSE to cause a rearrangement of the charge density 
that will reduce the electronic contribution to the electronic polarization when considering 
the polar phase. This can be understood as follows. The increase of the Jahn- Teller splitting 
goes hand in hand with the increase of the energy gap: the larger the gap, the smaller the 
dielectric constant is, i.e. the smaller the screening is. Now, let us consider the fixed 
ionic configuration of the paraelectric phase: the charge in the non-centrosymmetric spin 
arrangement can be thought as a "small" perturbation of the centrosymmetric one upon 
the application of the "internal" electric field. The electronic charge will respond to such 
a field, and each electronic state will change assuming a polarized configuration. If the 
gap is large, the "electric field" will hardly mix the electronic states in the valence band 
with the electronic states in the conduction band, since in second order perturbation theory 
approach the denominator will be of order of the band gap energy, so that the electrons 
don't polarize much, i.e. the charge distribution becomes more "rigid". In conclusion, we 
expect the electronic contribution to the electric polarization to decrease upon introduction 
of Fock exchange. This will be confirmed by our calculations. Note that the above reasoning 
is not appropriate for BFO where the polarization is mainly due to ionic displacements. 



B. AFM-E phase 

1. Structural properties 

Let us focus on the AFM-E phase, where the resulting symmetry is lowered by the 
spin configuration with respect to the AFM-A spin arrangements by removing the inversion 
symmetry. In Fig. [61 we show the relevant structural internal parameters, for the relaxed 
PBE and HSE structure, by considering the Mn-O-Mn-O-Mn chain (compare Fig. |6]and Fig. 



Hj). The Mn-0 short bond- lengths do not show significant differences between the parallel 
and antiparallel spin configuration in PBE as well in HSE. On the other hand, the long 
Mn-0 bond lengths are mostly affected: their difference, V - l ap in PBE is about 0.07 A 
and decreases to ~ 0.02 A upon introduction of exact exchange. At the same time, the 
angle changes: a p decreases while a ap remains almost equal to the PBE value. The results 
can be interpreted as follows: bare PBE is expected to overestimate hybridization effects 
between oxygen p-states and Mn d-states, therefore inducing a stronger rearrangement of 
ionic positions consistent with a "softer" structure when moving from, say a centrosymmetric 
A-type to a ferroelectric E-type phase. Viceversa, upon introduction of correlation effects, 
the reduced hybridization is expected to lead to a more "rigid" ionic arrangement. Indeed, 
this is evident when comparing the difference between a p - a ap , which drastically reduces 
upon introduction of HSE with respect to PBE. We recall that, ultimately, it is this difference 
that gives rise to the ionic polarization, as clearly shown in Fig. 2 b) of Ref. 49]. We can 
therefore anticipate that a reduction of the polarization will occur upon introduction of 
HSE, as discussed in detail below. What is worthwhile noting is that the Mn-Mn distances 
in PBE dramatically depends on their having parallel (d p Mn _ Mn = 3.98 A) or antiparallel 
spins (c?Mn-Mn = 3-87 A), so that d p Mn _ Mn — d ap In _ Mn = 0.11 A. However, this dependence 
is smoothed upon introduction of HSE, so that the difference in Mn-Mn distance strongly 
reduces to d p Mn _ Mn — d ap ln _ Mn = 0.03 A. In general, the marked (weak) dependence of the 
structural properties within PBE (HSE) is consistent with a strong (small) efficiency of the 
double-exchange mechanism, which ultimately relies on the p — d hybridization and hopping 
integral. 



2. Electronic and magnetic properties 

The band-structure of the AFM-E is quite similar to the A-type and is therefore not 
shown. However, there are some small differences which we comment on. As expected, the 
increase of the number of the AFM bonds of each Mn with its four neighbors associated with 
the change of the magnetic state going from AFM-A to AFM-E type results in a narrowing 
of all bands. This is further enhanced by HSE due to the reduced hopping upon introduction 
of exact exchange, as expected. Furthermore, the increase of the band gap is facilitated by 
the interplay of the crystal distortion, which is generally enhanced by HSE, with the AFM 



arrangement of spins. As expected, the energy gap is the largest in the AFM-E-HSE band 
structure, being now ~ 3 eV. The Ajt evaluated at S point, is also the largest in this 
case, being ~ 3.7 eV. Before turning our attention to the electronic polarization, we discuss 
the magnetic properties. First of all, we found that the AFM-E is more stable than the 
AFM-A by ~ 4 meV/cell in the HSE formalism. Note that this value has been obtained 
using the same simulation cell for both phases, therefore reducing the influence of numerical 
errors. Although the relative stability is still comparable with the numerical accuracy, it is 



indeed consistent with experiments. 
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In AFM-E-PBE, the Mn moment is 3.4 \i B which 
induces a small spin-polarization on the oxygen equal to ± 0.04 \lb- In AFM-E-HSE, the 
Mn moment slightly increases to 3.7 \i B and the oxygen moment slightly decreases to ± 0.01 
[Lb'- the increased localization of the Mn d states correlates with the increased Manganese 
spin moment and goes hand by hand with the decreased p-d hybridization and a decreased 
induced spin moment on oxygens. 



3. Ferroelectric "properties 

Finally, we calculated the ferroelectric polarization by considering the AFM-A as the 
reference paraelectric structure. The results show that the polarization (both the electronic 
and ionic terms) strongly reduces upon introduction of HSE. However, it is remarkable that 
the total P is still of the order of 2 fiC/cm 2 : this confirms HoMnC>3 as the magnetically- 
induced ferroelectric having the highest polarization predicted so far. Our estimate is in very 



good agreement with model Hamiltonian calculations. 89] The comparison between theory 



and experiments as far as the electric polarization is concerned is still a matter of debate 



Whereas earlier studies predicted negligible values for polycrystalline HoMnC>3 samples, 90] 



more recent studies for TmMnOs in the E-type (where the exchange-striction mechanism is 



exactly the same as in HoMnOs) point to a polarization which could exceed 1 [iC/cm 2 , 91| 
in excellent agreement with our predicted HSE value. 

The reasons why we expect a reduction upon introduction of HSE have been already 
discussed in previous paragraphs and can be traced back to the reduced p — d hybridization. 
As in the case of BFO, we disentangle the structural and electronic effects, by using the 
HSE (PBE) geometry with the PBE (HSE) functional (cfr. Table [TVT) . What we infer from 
these "ad- hoc" -built systems is that the use of HSE dramatically reduces the electronic con- 



tribution (cfr P e \ e in PBE and HSE(PBE)), i.e. reduced by ~ 2.5 fiC/cm 2 . Less important, 
though still appreciable, seem to be the ionic displacements: their dipole moment is reduced 
by ~ 1.5 fiC/cm 2 when comparing P tot in HSE and HSE(PBE). This is consistent with what 
shown approximately by P pcm . 



V. CONCLUSIONS 



In this work, we have revised the two workhorse materials of the exponentially growing 
field of multiferroics, namely BiFeOa for proper MFs and HoMnC>3 for improper MFs by 
using the screened hybrid functional (HSE). 

From our study, several important points emerge. For BFO: i) the structural, electronic 
and magnetic properties well agree with experiments; ii) the ferroelectric polarization agrees 
with reported values in the literature; iii) even if PBE allows the description of ferroelectric 
properties by opening an energy gap, it is by no means satisfactory in correctly describing all 
the properties on the same footing. On the other hand, HSE improves the PBE and LDA+t/ 
description; this is clearly shown by benchmark calculations using the most advanced and 
accurate state-of-the art GW+vertex corrections (which basically confirm the HSE results); 
iv) the previous comment, and very recent studies 40[] suggests that optical properties, so 
far not investigated at all by ab-initio calculations for BFO, can be properly addressed 
within HSE. v) finally, we note that the electronic polarization increases upon introduction 
of exact exchange. For HMO, we note that: i) the HSE results are in good agreement with 
experiments when available; ii) the Jahn- Teller effect is correctly described in agreement 
with experiment; iii) despite a reduction of the polarization value with respect to PBE, 
HoMnOa still shows the highest P predicted among magnetically-induced ferroelectrics. 

We have shown that introduction of "correlation" effects may both enhance the polar- 
ization or reduce it: the former effect will most likely occur for proper MFs, and the latter 
for improper MFs, e.g. magnetically driven. Note that an increase of HSE polarization 
with respect to LDA, for example, is also found by Wahl et aZ. (43 1 for BaTi03, a standard 
ferroelectric compound. Also for BiFe03, an increase of polarization using DFT+£7 has 



been noticed. 
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92] The increase of ferroelectric polarization when including a fraction 
of exact exchange and using the theoretical equilibrium volume has been reported also for 
simple ferrolectric compound such as KNb03.[ 



gquili 

•H 



One final comment is in order: although the HSE results certainly point towards a truly 
realistic description, it is still possible that, to some extent, the good performances of HSE 
may be material-dependent, i.e. the universal 1/4 fraction of the exact exchange may be 
not appropriate for some specific material. What is certainly true is that the predictive 
capability of HSE, combined with its nowadays affordable computational costs, make the 
functional an attractive choice for the study of a wide range of materials, from well-behaved 
insulators to doped semiconductors exhibiting magnetic ordering, multifunctional complex 
oxides of interest for many industrial applications therefore representing a very good starting 
point for materials design. 
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TABLE I: Dependence of structural paramenters, magnetic moment and electronic band gap with 



respect to PBE (this work), Bl-WC (Ref. 46j) and HSE (this work). a r h, a r h, V denote the 
lattice parameter, the angle, and the volume of BiFeG"3 in space group R3c, in the rhombohedral 
setting. Fe (2a) and O (66) Wyckoff positions are also reported. Bi is in the (2a) (0,0,0) position 
(not reported in the table). Long (1) and short (s) bond lengths of Bi-0 and Fe-0 are given. 
f3o-Fe-o is the O-Fe-0 obtuse angle (equal to 180° in the ideal octahedron), fipe is Fe magnetic 



moment. Experimental structural parameters are from Rcf. 
(measured on polycrystalline powder) is from Ref. 



651 ] , the experimental magnetic moment 



951 ]; the band gap is from Refs. 



66J471]. 
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HoMn0 3 (AFM-E) 
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5.18 (5.00) 7.41 (6.77) 


Mn 


6.82 (6.75) 12.10 (11.84) 





7.33 (7.41) 12.42 (12.61) 



TABLE II: Calculated charge and volumes according to Bader AIM partitioning for a fixed HSE 
geometry for BiFe0 3 and HoMn0 3 . Numbers in parenthesis are based on the HSE charge density. 
In all cases, we used the ferroelectric phase. 
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0.4617 


0.4616 


0.4622 


z 


0.6162 


0.6179 


0.6113 


2 8d(xyz) 








X 


0.3250 


0.3301 


0.3281 


V 


0.0550 


0.0573 


0.0534 


z 


0.1988 


0.2018 


0.2013 


Mn-0 (s-inpl) 


1.9271 


1.9022 


1.9044 


Mn-0 (1-inpl) 


2.2030 


2.2393 


2.2226 


Mn-0 (outpl) 


1.9518 


1.9546 


1.9435 


Q 


0.55 


0.61 


0.60 



a 19.47 19.71 18.77 

Mn-O-Mn (inpl) 143.849 142.839 144.081 
Mn-O-Mn (outpl) 141.060 140.586 142.462 



TABLE III: Structural parameters for HoMnC>3 optimized with AFM-A configuration in the Pnma 
unit cell. The lattice parameters are a=5. 8354, 6=7. 3606, c=5. 2572 A. Long and small Mn-0 dis- 
tances in c — a plane (inpl) and middle Mn-0 distance along the b axis (outpl). Mn-O-Mn angles 
in the c — a plane (inpl) and interplane Mn-O-Mn angles (with apical O along the b axis). 



BiFe03 Pjonic Pe/e Ptoi Ppem 

PBE 171.1 -65.5 105.6 87.8 

HSE 177.4 -67.1 110.3 103.4 

PBE(HSE) 177.4 -64.8 112.6 103.4 
HSE(PBE) 171.1 -67.9 103.2 87.8 
HoMn0 3 

PBE -0.6 -5.2 -5.8 -2.0 

HSE -0.3 -1.6 -1.9 -0.9 

PBE(HSE) -0.3 -3.1 -3.4 -0.9 
HSE(PBE) -0.6 -2.9 -3.5 -2.0 

TABLE IV: Ferroelectric polarization of BiFe03 and HoMnC>3, using PBE and HSE and Point 
Charge Model estimate is also given (P pcm ). Units are in fj,C/cm 2 . PBE (HSE) means a PBE 
calculation at fixed HSE geometry for the ferroelectric as well as for the paraelectric structure. 
Viceversa for HSE (PBE). For BiFeOs the polarization is along the [111] direction; for H0M11O3 
is along the c axis. The partial ionic and electronic contributions are also given. Obviously, the 
relative weigth of ionic and electronic contributions depends on the valence electronic configuration, 
e.g. Bi d in the core or in the valence. Therefore, given an electronic configuration, we are only 
interested in trends of Ptot and P pcm as far as the exchange-correlation functional is concerned. 




Paraelectric (R3c) Ferroelectric (R3c) 



FIG. 1: (Color on line) Paraelectric R3c and ferroelectric R3c phase of BiFe03 in the rhombohedral 
setting. Black spheres are oxygen atoms. The octahedron centered at one Fe atom is shown. Note 
the off-centering of the Fe atom and the corresponding octahedra distortion. 
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FIG. 2: (Color on line) PBE (left) and HSE (right) density of states for the ferroelectric phase 
of BiFeO"3. In black lines: (a) total DOS for one spin channel, (b) local Fe Dos for both spin 
components (minority spin are shown as negative), (c) local oxygens DOS for one spin component 
only. Additionally, in HSE (b), we show the Fe d DOS calculated using vertex corrections in W 
(see text) by red dashed line and the experimental spectral weight of d Fe taken from Ref.[3] by 
dotted blue lines; in (c), the magenta dashed line is Bi lone-pair DOS, i.e. Bi s states. The zero 
is set to valence band maximum. 




FIG. 3: (Color on line) 3D charge density difference (Ap = p — p ) isosurfaces for BiFe03. 
Yellow (grey) regions correspond to an excess of HSE (PBE) charge. Only the distorted octahedron 
at Fe site is shown. 



Paraelectric (AFM-A) 



Ferroelectric (AFM-E) 



FIG. 4: (Color on line) Upper part: perspective view of HoMnO"3; polyhedra surrounding the Mn 
atoms are also shown. Bottom part: paraelectric (left part) and ferroelectric (right part) spin 
configuration in the c — a plane. The electric polarization develops along the positive c axis. Red 
circles denote Ho atoms. 




FIG. 5: (Color on line) Band structure for the AFM-A phase of HoMn03 calculated along the 
symmetry lines of the orthorhombic Brillouin zone, for PBE (left) and HSE (right) using the 
relaxed structures. Red (green) dots refer to bands projected onto spin-up (-down) Mn atoms. 




Mn-Mn p =3.982 Mn-Mn p =3.944 

Mn-Mn ap =3.872 Mn-Mn ap =3.910 



PBE HSE 

FIG. 6: (Color on line) Schematic representation of the Mn-O-Mn-O-Mn chain in the c — a plane, 
corresponding to the up- down- down spin configuration. The structures are in scale and all struc- 
tural details are shown in the Figure (see text for further details). Distances are in A and angles 
are in decimal degrees (°). The left (right) part corresponds to the PBE (HSE) relaxed structure. 
The atomic displacements projected into the c — a plane which bring the PBE structure to the 
HSE one are shown by black arrows (left part). The spin configuration is shown only in the HSE 
structure by blue arrows. Blue (green) spheres are Mn (oxygen) atoms. Ho atoms are not shown 
for clarity. 



